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Abstract 

M, 

We study the chemical conditions at freeze-out associated with the production of 
strange-particles in Si-Au collisions at 14.6 GeV A. We obtain freeze-out chemical 
potentials and temperature, and determine the entropy as well as the final particle 
abundance. We also consider in detail the alternative evolution scenarios involving the 
hadronic gas and the deconfined phase. 
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1. Introduction. In the past years a number of experiments at 14.6 GeV A (AGS-BNL) 
and at 200 GeV A (SPS-CERN) have shown that the strange particle abundance differs 
Ph! considerably in nuclear interactions from that in nucleon-nucleon collisions |E|. Qualitative 
agreement with the experiment within the ARC particle cascade models has been reached at 
the AGS energy [§, ||, without need for any novel mechanisms of strange particle production 
enhancement, but at CERN energies such agreement has not been obtained [f| - moreover 
so far no microscopic model has emerged, capable to describe in all necessary detail the multi- 
strange anti-baryon abundance seen at CERN in the central rapidity region. Furthermore 
one finds significant difference in the strangeness flow comparing 14.6 GeV A with 200 GeV 
A collisions ||. 

The ARC results indicate transient presence of a clearly defined region of high density 
matter at central rapidity, seen in particular in the case of the Au-Au collision simulations. 
This suggests that a local equilibrium thermal model (fireball model) may also do well in 
describing these collisions and in this work we will assume simply that all central rapidity 
strange particles originate from a central, thermally equilibrated fireball — our analysis will 
allow in principle for chemical strange particle abundance deviations from full equilibrium 
here assumed. We implemented here several features which are more difficult to account for 
in the full microscopic simulation: 

• we incorporate the quantum gas nature for the hadronic particles; 

• we treat the entropy contained in the initial thermal state as a parameter and determine 
it as needed to describe the dat a; 
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• we explore the consistency of the description with the possible transient presence of a 
novel form of matter, such as is the locally deconfined quark-gluon plasma (QGP). 

The parameters of the thermal model are the temperature T and the chemical potentials 
\ii of the different quark flavors. In this discussion it is not necessary to distinguish the u 
and d flavors and we therefore consider the light quark chemical potential /i q and the related 
fugacity A q = exp(/i q /T). Since three quarks together have baryon number one, the baryon 
number fugacity is Ab = A q , that is /xb = 3/x q . A finite value of /is indicates that we have a 
finite baryon density. 

We will study in some detail strange particles which contain strange quarks (or anti- 
quarks), associated with strange quark fugacity and chemical potential: A s = exp(/x s /T). At 
relative particle abundance equilibrium the value of fi s is fixed within our approach by the 
requirement that strangeness is conserved. In any hadronic-gas (HG) picture of the reaction 
leading to high baryon density, there exists typically a considerable asymmetry in the strange 
and anti-strange particle phase space and hence the requirement of strangeness conservation 
normally leads to a considerable value of /i^ G > 0. An exception from this rule exists and 
is discussed in detail elsewhere || |7J. On the other hand, since in the deconfined QGP 
state both s and s quarks are free, one always has fJ$ = 0, assuming that strangeness is 
conserved, < s — s >= 0, on the time scale of A-A collisions. Said differently, observation of 
/i s 7^ (that is strange-quark fugacity A s ^ 1) strongly suggests that these particles emerge 
from a final state which is a confined HG phase. 

In the analysis of the collisions of S-ions at 200 GeV A with different nuclear targets 
carried out at CERN-SPS one finds in the framework of the thermal model || |§ that the 
strange-quark fugacity is A s ~ 1 (or said alternatively, the strange-quark chemical potential 
/4 ~ 0). Furthermore, considering particle production one finds that these collisions are lead- 
ing to very entropy rich state fllOfl . Finally, also the strange-particle phase space occupancy 
is close to equilibrium, (7 S ~ 1) || — this suggests the presence of an unusually effective 
strangeness production mechanisms expected only in the QGP phase. These three results 
strongly suggest that QGP phase is formed in collisions at 200 GeV A, but considerable 
additional effort is needed to show this to be the case. Importantly, /i s = is found, despite 
the expectation that in (a slow) hadronization process of the QGP the memory about the 
transient, deconfined state would be erased and hadronic particle abundance would reflect 
conditions of an equilibrated HG. Thus the value /i s = also suggests (if we presume forma- 
tion of QGP like state) that at CERN the hadronization is sudden and that therefore certain 
hadronic observables could carry considerable information about the primordial phase. One 
of our principal results here is a demonstration that as the energy decreases from ^/snn = 8.8 
GeV, at CERN, to a/snn = 2.5 GeV, at AGS, this situation changes with the freeze-out well 
described by equilibrium HG. This may be considered to indicate that at AGS: 

• the QGP phase was not formed at all, or, 

• that complete re-equilibration occurs when the primordial QGP hadronizes. 

While the freeze-out chemical potentials at CERN energies are well known, the freeze-out 
temperature remains yet to be determined, despite the relatively large sample of data: the 
results indicate that the hadronization of the primordial phase is sudden, and one has to 
permit non-equilibrium abundances of diverse hadronic particles. This introduces additional 
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off-equilibrium parameters which need to be determined along with the freeze-out tempera- 
ture . On the other hand there seems to be no need to introduce non-equilibrium particle 
abundance parameters (other than strange phase space occupancy factor 7 S ) at AGS energies 
and hence we can determine in this work the freeze-out temperature applicable to AGS - 
we note that the inverse slope parameters of particle spectra, interpreted within the fireball 
model are the blue shifted freeze-out temperatures Tf, the Doppler shift arising from the 
collective flow phenomena |TT] . 

In this work we will employ state of the art equations of state for both HG and QGP 
fireballs: 

• for HG we use a refined Hagedorn approach which accounts for interactions by treating 
all hadronic resonances as independent particles. We include numerically the known 
resonances, considered usually to be Boltzmann particles f6j, but in the present work 
which reaches to lower temperatures T < m n , where m n is the pion mass, pions and all 
mesons are treated as bosons and we also treat nucleons and A-resonances as fermions. 

• for QGP we take the analytical form of an ideal quark-gluon gas, but amend the number 
of degrees of freedom to reflect the reduction due to perturbative QCD-interactions, 
and we introduce non-perturbative thermal mass for quarks and gluons [12[ . 



2. AGS data. We will discuss here Si-Au results obtained at 14.6 GeV A by the experiment 
E802/E859 at BNL, which were presented at the QM'93 meeting [ 1I3L O] , in particular: 

• the ratio K~/K + = 0.22 ±0.02, which is constant in the central rapidity range 1 < y < 
1.6, consistent with the hypothesis that the dominant particle source is a fireball. This 
ratio varies strongly in target and projectile rapidity regions, as can be expected due to 
secondary production mechanisms and the presence of baryon-rich reaction fragments; 

• the central rapidity ratio A/A = (2.0 ± 0.8)10 -3 averaged over the acceptance of the 
spectrometer (1.1 < y < 1.7). 

We will now use these strange particle ratios to determine the values of chemical potentials 
at freeze-out; we have R [7], |8|, M: 



A 
A 



K 2 K' (!) 



K 

K+ 



T7T = W , (2) 



We obtain 0: 



A _ (EVKY^ _ f +0.20 

s ~ TaTaT 7 ^ i -°- 15 ' ( ) 

(K+/K-) 1/6 / +0.39 

= (a/a) 1/6 = I -°- 26 ' 
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These two results are considerably different from the result obtained at CERN energies: 
A s = 1.03 ±0.05, A q = 1.49 ±0.05. The AGS result for A s differs by more than four standard 
deviations from the CERN value — we note that the high-power roots entering above analysis 
reduce the significance of the large errors in the data. We further note the relatively high 
value of A q which is more than twice as large as that of the CERN experiments and suggestive 
of considerable baryon density reached at AGS conditions. 

An interesting implication and to some extent a verification of the thermal fireball model 
is that the two particle ratios we consider directly determine the p/p ratio — note that 
Eqs. (D, © imply: 



^ = A- 6 = ^| T = (0.44±0.2)10- 3 . (5) 
We note fl3| that at central rapidity y — 1.3 the p-yield is given to be dN/dy ~ 10~ 2 , which, 



when combined with the central proton yield dN/dy ~ 20, leads to p/p = (0.5 ± 0.2) 10 3 
(where we estimate the error from the results shown in |13|]), which agrees in a remarkable 



way with the result we obtained above. Since within the error bar the ratio varies by factor 
two, it is also consistent with the other experimental results reported [[15[]. We stress that 



this agreement of the fireball model with experiment may be fortitouse, since the p data are 
still preliminary, while the p-yield is strongly influenced by the target spectator tail. We will 
therefore return to a full discussion of the validity of the fireball model when the Au-Au 
results become available. 

3. Theoretical interpretation. A natural question arising here is in how far these values 
of A s , A q are consistent with the expected properties of the hadronic state and what they 
tell us about other freeze-out conditions. First, let us note a powerful constraints for our 
analysis: the three thermal parameters T, A s , A q are constrained by the requirement that the 
net strangeness of the source be zero — roughly speaking, the sum of hyperons and anti-kaons 
must equal the number of kaons, (we of course account for all other strange particles) and the 
values found in Eqs. (|3]), (^) are consistent with this constraint. When implementing in detail 
the strangeness conservation constraint we allow for one off-equilibrium parameter: since 
strange quark flavor is produced relatively slowly in hadronic interactions (but strangeness 
produced is easily redistributed among different degrees of freedom), it is customary || 
to weight with a factor 7" all hadrons containing n strange or anti-strange quarks. This 
strangeness phase space occupancy factor enters in some theoretical results 0,thus we will 
allow in the discussion below a range of values: on the low side we take 7 S = 0.2, which is 
the value typical of p-p collisions, and the upper limit will be 7 S = 0.7, which is found in 200 
GeV A experiments. The equilibrium value is 7^ = 1. 

A parameter of all thermal models is the specific entropy (entropy per baryon) content 
of the source: once the initial thermalized state is established, its evolution to the freeze-out 
condition is proceeding nearly without entropy production as we have found studying the 
approach to chemical equilibrium for (isolated) fireballs fl6|| . A most remarkable property 



of the HG (and also QGP) equations of state of hadronic matter which we use here is that 
at a given value of specific entropy per baryon (given strangeness conservation constraint) 
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corresponds to a nearly constant ratio: 

Rx = A q /A s . (6) 

In the QGP phase this is well understood [§, Q - a dimensionless quantity, such as is 
S/B can only be a function of dimensionless objects, and since A s = 1, it is only a function 
of A q . In HG this finding which implies a large degree of cancellation, was for us totally 
unexpected. In Fig. [I] we show for each value of specific entropy 10 < S/B < 30 the range 
of R\ = A q /A s which is possible when we vary temperature in the wide range T G (100, 210) 
MeV and also allow for a wide range in strange particle occupancy 7 S G (0.2, 0.7) (lower 
bars are for 7 S = 0.2 and upper ones for 0.7). We have interpreted in Eq.(|2]) the kaon 
ratio in terms of fugacities. This implies R\ = 2.12 ± 0.12 and we show (right two vertical 
lines in Fig. |I|) that this corresponds to the specific entropy at freeze-out S/B = 13 ± 1 (for 
7 S ~ 0.4). This implies after all resonances have decayed a final state ratio of mesons to 
baryons R m /b — 1.6 ± 0.2 , which is probably compatible with the somewhat lower value 1.2 
reported 0; since in our calculus we account for the enhancement of low p± pions, arising 
from resonance disintegration. In Fig. [I] we also show the value of R\ = A q /A s ~ 1.5 
associated with the 200 GeV A results (left two vertical lines in Fig.|l|) — we see that the 
HG type interpretation at CERN would imply a specific entropy content of 24 ± 5, while 
the high particle flow seen indeed points towards S/B = 40 more appropriate for a rapidly 
disintegrating QGP phase [[H| • 



In order to determine the freeze-out temperature Tf, we consider (see Fig. ^) the ratio 
A/A at fixed S/B as function of T. Again we find that there is residual 7 S dependence — the 
thick lines are for the upper limit 7 S = 0.7, the thin lines for the lower limit 7 S = 0.2. The 
short dashed (upper) lines are for S/B = 14, the long dashed (lower) ones for S/B = 12 and 
the solid lines are for S/B = 13. The experimental range A/A G (1.2, 2.8) 10 -3 thus leads to 
the range Tf = 127 ± 8 MeV, and we indicate the experimental and theoretical ranges of the 
values by the error cross in Fig. |[ Given the thermal freeze-out conditions we can obtain 
the freeze-out baryon density pf = (0.7 ± 0.3)ptv, where pn = 0.16 fm~ 3 . 

We now turn to the two questions: 

1) what are the initial thermal conditions when the fireball is formed? 

2) what differences, if any, are there between HG and QGP phases? 

In Fig. [5] we show in the pe, T plane for the HG case (thick lines) and QGP case (thin 
lines) the different possible histories of the collision. The dashed lines correspond to a fixed 
energy per baryon 2.55 GeV A in the CM frame, while the solid lines are for fixed specific 
entropy per baryon S/B = 13. Where the solid and dashed lines meet (within a trapezoidal 
region determined by one unit error in entropy and an error of 0.15 GeV in CM energy) we 
have a consistency between initial energy and final state entropy, thus presumably these are 
the initial values of thermal parameters. The initial HG state (for S/B = 13, E/B = 2.55 
GeV) has an energy density e^ G = 2.3 GeV fm~ 3 , and baryon density p^ G = 6pN- The 
initial pressure is P^ G = 0.3 GeV fm -3 , which is normally considered a value typical of 
the deconfined phase. In the QGP phase we find accordingly: £^ GP = 1.2 GeV fm -3 , 
P?° P = 2.9pAr, P^ GP =0.39 GeV fm~ 3 . Surprizingly, the QGP is the more dilute phase at 
these conditions, but it can be argued that the greater stiffness of QGP, as compared to HG, 
reduced the nuclear compression reached in collision. 
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Remarkably, we find a very large difference in /xb which takes an initial value /i^ G = 
440 ± 40 MeV in the HG scenario and u$ GP = 910 ± 150 MeV in the QGP case — not 
shown in the Fig. |3] is that the strangeness conservation requirement leads in the HG to an 
initial value /i^ G = 0, just as is in the case of the QGP. We see in Fig. |3| that the initial 
temperatures T for QGP and HG scenarios are practically equal. For E/B = 2.55 ± 0.15 
GeV and S/B = 13 ± 1 we have T = 190 ±30 MeV. We recall that the effective temperature 
seen in the m± particle spectra is nearly To || 11]. Thus presence of high inverse slopes 



in high m± spectra with T ~ 225 MeV [13|] is consistent with our considerations and both 



equations of state used (HG and QGP). Perhaps a small edge can be given to the QGP case, 
since one should allow in principle for some entropy production in the transition to HG and 
for S/B = 11 one would indeed reach with the central value to T^ GP ~ 215 MeV. 

The QGP fireball at S/B = 13 (thin solid line) evolves practically at fixed A q = 4.8 
and A s = 1. For the HG fireball at fixed S/B = 13 (thick solid line) there is a strong 
variation in both these fugacities but as noted in Fig. ^] the ratio R^ G (= A q /A s ) = 2.17 
remains practically constant. The 'experimental' cross is set at /iB,f = 485 ± 70 MeV and 
Tf = 127±8 MeV corresponding to our determination of freeze-out conditions (with /i s ,f = 68 
MeV). This appears near the region where the two evolution curves meet (QGP and HG) 
but in fact both miss each other greatly in a three dimensional display, since the ratio 
R^ GP = A q /A s = 4.8 on the QGP path-line and = 2.17 on the HG path-line. 

We have also studied where a possible phase transformation from QGP to HG may be 
if QGP is produced. Since we are in relative low pressure domain, the pressure of the QGP 
phase is impacted by the vacuum bag constant B contribution. We took B — 0.1 GeVfm -3 
(corresponding to B 1 ^ = 170 MeV). From the resulting phase diagram we obtained that 
in the event that a first order phase transition is occurring, it would be at P = 0.04 ± 
0.01 GeVfm -3 and this determines points on the two evolution paths shown, connected by 
nearly vertical lines — we note that there is minor re-heating occurring while the baryo- 
chemical potential drops by 15% — however, the major re-equilibration is in the jump from 
A s = 1 in the plasma to A s ~ 1.7 in the HG phase (as discussed above this is still amplified 
by comparing the greatly different values of R\ in both phases). Thus the 'short' connection 
between the QGP to HG paths would be considerably stretched in full three dimensional 
display, reflecting on the need to well re-equilibrate the matter in transition. Therefore we 
believe that there is not a slightest remembrance of the possible QGP history in the AGS 
data considered here. 

However, it is likely that the presence of a primordial phase can be in principle seen in 
hadronic particle spectra due to the great difference in the R\ ratio between QGP and HG 
cases — the difference is by factor 2.2. Thus e.g. the ratio of K + to tt + particle abundance 
at fixed m± should increase at large m± supposing that the high m± particles are directly 
emitted from the possibly formed QGP phase. Since in the fireball model 



K 

7? 



(7) 



we should see that the high m± ratio decreases just by factor 2.2 as the high m± QGP 
primordial contribution are overcome by hadronization yields, and at still smaller m± by 
even more as the pion abundance is diluted by the A decays (an effect which we can fully 
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account for given the model for the freeze-out conditions we developed). 



4. Conclusions and outlook. We have reached a full and satisfactory description of the 
freeze-out conditions at AGS energies, employing thermal fireball model. We have alluded 
to a number of consistency checks, e.g. have mentioned that the meson to baryon ratio in 
the final state R m /b = 1.6 ±0.2, and also noted that the ratio p/p = 0.5 10~ 3 is of the correct 
magnitude. We note that given the freeze-out conditions we can compute the abundances 
of most hadronic particles in the final state, and we intend to develop a more detailed and 
comprehensive analysis once more precise input data become available. 

Should QGP be formed at these low energies, we hope that akin to the situation at 
200 GeV A, one will be able to see considerable anomalies in abundances of multi-strange 
anti-baryons. These are produced extremely rarely in the HG environment and normally one 
would not believe that they can be observed for beam energies at about 14.6 GeV A. However, 
should at AGS indeed the QGP phase be formed, their abundance would be only little 
reduced as compared to the 200 GeV A situation, while the experimental sensitivity would 
be greatly amplified, since the entropy content and hence the final state multiplicity is much 
smaller at these lower energies. Given the theoretical results we have presented regarding 
the initial thermal conditions reached in the central fireball in 14.6 GeV A collisions, it is 
hard not to envision QGP formation at AGS and hence we would like to encourage strongly 
that a search be made for these hadronic observables [[H]] of the QGP phase. Thus we are 
in particular interested in central rapidity data about (multi-) strange anti-baryons which 
we hope could help resolve the question if there has been formation of transient deconfined 
phase. In particular we may be able to determine the value of 7 S ; we need to look at the 
ratios: 



7s 



%R> 



1 A + S° 

2 p 



A + S c 



1 n- 

2 1= 



Resonance decays influence these yields, and hence the above relation holds approximately, 
but increasingly exactly as the strangeness content increases. If QGP is formed and the 
strangeness phase space were fully saturated, we would expect that the R§ ratio is indeed 
greater than unity. Due to the large value of A q we found, if QGP is formed with 7 S ~ 0.7 
and R\ = A q /A s = 4.8 governing emission at high m±, we obtain R§ ~ 3. Using the HG 
values 7 S ~ 0.15, R\ = 2.1 one finds a yield typical for p-p reactions, with R$ ~ 0.3. We 
recall here that for the CERN experiment the corresponding discussion addressed the ratio 
of yields of S/A and the high value of 7 S ~ 0.7 found there arises from the remarkable result 
that R- s is 0(1) at 200 GeV. 

Even though we have here discussed only the Si beam experiments at 14.6 GeV A in- 
volving collisions with the heaviest target (Au) as is appropriate for the here considered 
fireball model, our analysis can be easily continued to the case of Au-Au collisions at 11-12 
GeV A [D, [18| . Because of the projectile-target symmetry the CM specific energy (assuming 
full stopping or equal baryon and energy stopping) is for the generally considered value of 
11.7 GeV A exactly the value we assumed here, namely v / snn = 2.55 GeV. We note that a 
slight increase in the ratio of K + /K~ from 4.5 to perhaps 5 would imply a slight drop in the 
specific entropy (see Eq. (0) and Fig. [I]), perhaps by 1/2 unit. This causes only incremental 
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changes in all freeze-out results and thus there will be little change in final state particle 
abundances. However, we find, within the framework of HG model, a considerable sensitivity 
of the initial conditions (baryon density, pressure and energy density), seen as function of 
the entropy content: a drop by one unit of S/B from 13 to 12 implies doubling of all these 
initial values. Thus it is well possible that much more extreme conditions are reached in 
the Au-Au collisions as compared to Si-Au, even if the reaction products and multiplicity 
appear to be similar for both systems. 

In summary: we have used that results on strange particle ratios A/A and K~/K + from 
AGS in a fireball model framework and have determined freeze-out values of temperature 
(T f = 127 ± 8 MeV), and chemical potentials (/i B ,f = 485 ± 70 MeV and /i Sjf = 68(+20, -16) 
MeV, final state specific entropy S/B — 13 ± 1 accompanied by meson to baryon particle 
ratio i? m /b = 1-6 ± 0.2, and also within the HG and QGP models, the (extreme) initial 
conditions prevailing when the fireball is formed, in particular To = 190 ± 30 MeV. 
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FIGURE CAPTIONS 

Figure 1: For a given specific entropy per baryon S/B we show the range of allowable 
values A q /A s values in HG phase, when temperature is varied between 100 and 210 MeV and 
the strangeness occupancy factor 7 S between is 0.2 (lower values) and 0.7 (upper values). 
Strangeness conservation constraint is imposed. See text for discussion of vertical lines. 

Figure 2: The ratio A/A as function of temperature, at fixed specific entropy and for 
7 S = 0.7 (thick lines) and 0.2 (thin lines). Solid lines S/B = 13, short dashed S/B = 14, 
long dashed S/B = 12. 

Figure 3: f^B^T plane with thick lines: HG model; thin lines: QGP model. Dashed 
lines: fixed energy per baryon E/B = 2.55 GeV, upper for HG, lower for QGP. Solid lines: 
fixed specific entropy S/B = 13, upper for QGP, lower for HG model. Trapezoidal regions 
enclose the initial condition (see text). Solid line connecting QGP and HG with dotted lines 
left/right: possible phase transition at P = 0.04 ± 0.01 GeV/fm 3 . Note that continuations 
of S/B = 13 lines beyond transition or hadronization points are dotted. 
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FIGURE 3 



